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Abstract 

Grand rotation curves (GRC) within ~ 400 kpc of M31 and the Milky Way were constructed by 
combining disk rotation velocities and radial velocities of satellite galaxies and globular clusters. The GRC 
for the Milky Way was revised using the most recent Solar rotation velocity. The derived GRCs were 
deconvolved into a de Vaucouleurs bulge, exponential disk, and a dark halo with the Navarro-Frenk-White 
(NFW) density profile by the least fitting. Comparison of the best-fit parameters revealed similarity of 
the disks and bulges of the two galaxies, whereas the dark halo mass of M31 was found to be twice the 
Galaxy’s. We show that the NFW model may be a realistic approximation to the observed dark halos in 
these two giant spirals. 

Key words: galaxies: dark matter — galaxies: individual (M31) — galaxies: The Galaxy — galaxies: 
rotation curve 


1. Introduction 

Rotation curve of the Andromeda galaxy M31 (NGC 
224) has been obtained in details and are used for deter¬ 
mining the mass distribution in the disk and dark halo 
(Sofue et al. 1999; Carignan et al. 2006; Chemin et 
al. 2009; Corbelli et al. 2011). The dark halo mass be¬ 
yond the disk has been derived using kinematics of satel¬ 
lite galaxies (e.g., Metz et al. 2007; van der Marel et al. 
2008; Tollerud et al. 2012). Globular clusters are also 
used extensively for the mass determination of dark halo 
(Veljanovski et al. 2010) as well as for inner kinematics 
(Galletti et al. 2004). 

In our recent work on the Milky Way Galaxy, we con¬ 
structed a running averaged rotation curve for a wide re¬ 
gion from the Galactic Center to outermost dark halo, 
which we called the grand (or pseudo) rotation curve 
(hereafter, GRC: Sofue 2012, 2013). In the present paper, 
we revise the GRC by adopting the most recent value of 
the Solar rotation velocity, Vq = 238 km s“^, from VERA 
observations (Honma et al. 2012). 

We apply a common method to construct GRCs to M31 
and the Milky Way by combining disk rotation curves and 
radial velocities of satellite galaxies and outer globular 
clusters. We deconvolve the GRCs into bulge, disk and 
dark halo, and determine the dynamical parameters. We 
compare the results for M31 and the Milky Way, and dis¬ 
cuss the similarity of the two galaxies. It is emphasized 
that the NFW model indeed works in largely extended 
dark halos of real galaxies. 

2. Grand Rotation Curves 
2.1. MSI 

The rotation curve of the disk of M31 has been obtained 
by many authors using HI, CO and Ha line kinematical 


data as listed in table 1. For the disk region, we simply 
adopt the averaged rotation velocity V {R) as a function 
of radius R from these data. For the halo region beyond 
the disk, we employ kinematics of non-coplanar objects 
orbiting M31, which include satellite galaxies and globular 
clusters (table 1). We adopt a distance to M31 of 770 kpc, 
and the systemic velocity of 300 km s“^ (Courteau and van 
den Bergh 1999) 

We assume that the distribution of non-coplanar ob¬ 
jects is spherical around M31 and the velocity vectors are 
random. We define a pseudo rotation velocity, V, as the 
velocity that yields the enclosed mass M {R) within a ra¬ 
dius R from the galaxy’s center by 

M{R) = V^R/G, (1) 

where G is the gravitational constant. The pseudo rota¬ 
tion velocity is here replaced by the Virial velocity for an 
ensamble of particles orbiting around the center of mass 
in random orbits (Limber and Mathews 1960; Bahcall and 
Tremaine 1981). If the orbits are random, the velocity is 
evaluated by 

V = ^3 < >, (2) 

where Vz is line-of-sight velocity of each object (Limber 
and Mathews 1960), and < > is the mean of squares 

of Vz in the ensamble. The factor was multiplied for 
correction for the degree of freedom of random motion 
(Limber and Mathews 1960). 

If the orbits are not random, the above estimation may 
be corrected for the orbits shapes by 

V = G^/?><vl >. (3) 

Here, the coefficient C is a factor ranging from 1% ~ 
0.7 for radial orbits to 37r/8 ~ 1.2 for circular orbits 
(Bahcall and Tremaine 1981), depending on the shapes 
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Table 1. References to the data in figure 1 


Rectangles at A < 10 kpc 

Grey circles at i? < 32 kpc 

Grey open circles linked by line 
Black-grey circles linked by line 

Grey reverse triangles linked by line 

Disk RC (CO) 

ibid (combi: HI, CO, opt) 

ibid (HI) 

ibid (HI, CO) 

ibid (HI) 

Loinard et al. (1995) 

Sofue et al. (1981, 1999) 
Carignan et al. (2006) 
Chemin et al. (2009) 

Corbelli et al. (2010) 

Rectangles at R > 40 kpc 

Triangles with bars 

Reverse triangle at R > 40 kpc 

Open circles with bars 

galaxies around M31 

ibid 

ibid 

Globular clusters 

Metz et al. (2007) 
van der Marel et al. (2008) 
Tollerud et al. (2012) 
Veljanovski et al. (2014) 


of orbits as well as on the assumed potential and distri¬ 
bution function of objects (Evans et al. 2003). 

Since, it is rather the potential that we want to deter¬ 
mine in this paper, and the other factors are unknown, we 
here assume that C = 1. It should be mentioned that the 
obtained result in this paper depends on this correction 
factor in the sense that the dark halo mass is proporional 
to Namely, the mass could be changed by a factor 
of ^ 2. Hence, the given errors in this paper are statis¬ 
tical values, not considering the systematic effects by the 
assumption. 

For objects whose three dimensional coordinates are un¬ 
known, the distance from M31’s center R is related to the 
projected distance rproj by 

I? = Fproj • (4) 

We adopt thus calculated pseudo rotation velocity as the 
rotation velocity. 

If there exists systematic rotation, as observed for the 
globular clusters around MSI (Veljanovski et al 2014), a 
correction is necessary as described in Appendix. We as¬ 
sume that the rotation axis is parallel to the disk’s rotation 
axis, i = 77°, so that the correction is small and equation 
3 holds also for the globular cluster system. 

Figure 1 shows compiled rotation velocities for the disk, 
and pseudo rotation velocities of satellites and globular 
clusters within radius 500 kpc from M31. Using these ve¬ 
locities, we calculate running averaged velocities in each 
bin at logarithmic interval between R and 1.2R with R 
starting from i? = 0.1 kpc by applying the Gaussian run¬ 
ning averaging procedure described in Sofue (2012, 2013). 
Note that we here use the root of the mean, ^y<v^> as 
the pseudo rotation velocity, because the squared velocity 
is more directly related to the Virial mass according to 
equation 3, whereas we took < Vz > in the ealier work. 

The pseudo rotation velocities are plotted in figure 2, 
where long bars represent modified standard deviations 
among the used data in each bin, and short bars are 
modified standard errors. Here, the modified deviation 
and error are defined by Sd = SV ^ I2V = Sd:v‘^l‘2V 
and Se = Sf.,yil2V, recalling the propagation of the de¬ 
viation (error) of and V by the derivative relation 

SV^ = 2VSV. Here, sa-.v^ = (< (3v^ - and 

Se:y 2 = are the standard deviation and error 

around in each radius bin, respectively, and N is the 




Fig. 1. Rotation velocities of the disk and pseudo rotation 
velocities of non-coplanar objects in M31 (a) in linear and (b) 
semi-logarithmic scaling. References to the data are listed in 
table 1. 

number of data points in the bin. 

These definitions were employed for the reason why we 
chose the square value, uf, as the independent variable, 
instead of the linear value, Vz- We call thus obtained plot 
of V the grand rotation curve (GRC). 

2.2. The Milky Way 

We revise the GRC of the Milky Way obtained in our 
earlier work (Sofue 2013) by adopting the most recently 
determined value of the Solar rotation velocity, Vo = 238 
km s“^at Rq = 8 kpc, from VERA observations (Honma 
et al. 2012) in place of 200 km s“^in the earlier paper. We 
also adopt t he same correction factor for velocity as for 
M31, V = y/3 < vl >, in place of <Vz>. In the earlier 
work, the degree of freedom of motion was assumed to be 
two, considering that each particle has transverse velocity 
same as the radial velocity. 

The revised GRC is shown in figure 3, where used data 
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Fig. 2. GRC of M31, showing running averaged values of 
pseudo rotation velocities in (a) linear and (b) logarithmic 
scalings. Long and short bars represent modified standard 
deviations and errors, respectively. 




Fig. 3. New GRC of the Milky Way in (a) linear and (b) 
semi-logarithmic scalings. Grey dots are the data for indi¬ 
vidual objects, and grey line with triangles shows our earlier 
GRC (Sofue 2012, 2013). 


Fig. 4. New GRC of the Milky Way (full line) compared with 

that of M31 (grey line) in (a) logarithmic and (b) linear scal¬ 
ings. The bars are modified standard errors. 

are superposed by grey dots, among which a few objects 
with Vz ~ 400 km s“^have been removed from the analysis. 
The new GRC exhibits higher velocities than our eariler 
result shown by the grey line with triangles. The higher 
disk velocities are due to the adopted larger value of Vb = 
238 km s“^as well as to higher velocities for non-coplanar 
objects because of here adopted correction of a/S instead of 
•\/2. Also, the presently employed average of the squared 
velocity uf lead to slightly higher mean velocity compared 
to the linear mean of Vz in the previous work. 

In figure 4 we compare the new GRC of the Milky Way 
with M31. The two GRCs show a remarkable similarity 
in the halo regions, indicating similarity dark matter dis¬ 
tributions. 

2.3. Characteristics of the GRC 

The general characteristics of the obtained GRCs for 
M31 and the Galaxy may be summarized as follows: 

(i) Rotation curves in the bulges look typical in the 
plot up to 40 kpc, but the enlarged curves in logarith¬ 
mic presentation exhibits that they have inner structures. 
Particularly, the innermost curve of the Milky Way was 
shown to be composed of multiple bulges with steep con¬ 
centrations than the de Vaucouleurs law (Sofue 2013). 

(ii) The disk rotation curves are nearly flat, showing 
broad maxima at i? ~ 10 to 20 kpc in both galaxies. 

(iii) Beyond i? ~ 30 kpc the rotation velocity declines 
smoothly until the edge of the dark halo. The slope is, 
however, milder than the Keplerian law, implying that the 
halo cannot be represented by the Plummer or exponential 
type potentials that require steeper mass concentration. It 
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will be shown in the next section that the NFW model is 
a good approximation to represent the observed GRCs in 
the halos. 

3. Deconvolution of GRC into Bulge, Disk and 
Dark Halo 


We assume that the GRCs shown in figures 2 and 3 are 
composed of bulge, disk, and dark halo contributions as 

V{Rf = VURf + Vd{Rf + V^{Rf, (5) 

where V{R), Vh{R), 14(i?), and Vh{R) are the rotation 
velocity at galacto-centric distance R, and those for the 
bulge, disk and dark halo, respectively. By fitting the 
GRC with a model rotation curve using the the least 
method, following the method presented by Sofue (2012, 
2013), we searchd for the best-fit parameters of the bulge, 
disk and dark halo. 


3.1. Bulge 

The bulge is assumed to have the de Vaucouleurs (1958) 
profile for the surface mass density as 

,1/4 


Eb(r)=Ebeexp -K^{r / 


( 6 ) 


where k, = 7.6695, Ebe is the surface mass density at the 
half-mass scale radius R = ah. The total mass is calculated 
by 


poo 

Mb = 27r/ rEb(r)dr-= lyObEbe, 
Jo 


( 7 ) 


with r] = 22.665 being a dimensionless constant. The cir¬ 
cular rotation velocity is then given by 


Rb(i?) = VGM^{R)/R. 


( 8 ) 


In the fitting procedure, Mb and Ob are taken as the 
two free parameters. The bulge of our Galaxy was shown 
to be composed of multiple bulges with exponential den¬ 
sity profiles, whereas the de Vaucouleurs law rather fails 
to reproduce the innermost rotation curve (Sofue 2013). 
Hence, the present analysis will be not accurate enough 
for the discussion of the bulge in the Milky Way. 


3.2. Disk 


The galactic disk is approximated by an exponential 
disk, whose surface mass density is expressed as 

Sd(i?) = Soexp(-i?/ad), (9) 

where Eq is the central value and ad is the scale radius. 
The total mass of the exponential disk is given by 

nOO 

Md = / 27rrEddr = 27rEoab. (10) 

Jo 

The rotation curve for a thin exponential disk is expressed 
by 

Vd{R) = ^GMd/ad^{X), (11) 

where X = R/ad, and ^{X) is the expression obtained by 
Freeman (1970) for a flat exponential disk. As the two 
free parameters we chose Md and Od. 


3.3. Dark Halo 

For the dark halo, three mass models have been so far 
proposed: the semi-isothermal (Begeman et al. 1991), 
NFW (Navarro, Frenk and White 1996), and Burkert 
(1996) models. The outermost rotation curves in figures 

5 and 6 are not flat at all, so that the isothermal model is 
not a good approximation. Since the NFW and Burkert 
models are essentially the same except for the very central 
part, we here adopt the NFW profile. The NFW density 
profile is expressed as 

p{R) = Po/[X{l + X)% (12) 

where X = R/h, and po and h are the representative 
(scale) density and scale radius of the dark halo, respec¬ 
tively. In the fitting procedure, we chose po and h as the 
two free parameters. 

The enclosed mass within radius R is given by 

Mh(R) = dTTpoh^ {ln(I + X)-X/{1+ A)} . (13) 

The circular rotation velocity is given by 

Vd{R) = ^GMd{R)/R. (14) 

3 . 4 . Fitting result by the least method 

Applying the fitting method described in Sofue (2012) 
to the GRC of M31 and the Galaxy, we searched for the 
best-fit parameters of the bulge, disk and dark halo. The 
free parameters are Mb, Ob, Md, Od, po and h. Fitting 
radii were taken to be i?i = 0 to R 2 = 20 kpc for the bulge, 
0 to 40 kpc for disk, and 1 to 385 kpc for the dark halo. 
The outer boundary for the halo fitting corresponds to 
the half distance between the two galaxies. Figures 5 and 

6 show the thus obtained fitting results compared with 
the used GRCs. Table 2 shows the best-fit parameters for 
individual mass components. 

Figures 7 and 8 show the behaviors of around 

the least values, where x is defined by x^ = S[14(Ri)^ ~ 
Vo{Ri)‘^]/s‘^ with o and c abbreviating observed and calcu¬ 
lated values, respectively, and N is the number of fitting 
points. For the dark halo we present the total dark mass, 
Mb : 385 , corresponding to h and po. Since the fitting areas 
and N for bulge, disk, and halo are different, the minimum 
X values are different among the components. The error 
of each fitted parameter was evaluated as the range that 
allows for an increase of the x^ value by 10% above the 
least value. 

4. Discussion 

4 . 1 . Summary 

We constructed GRCs of the Andromeda galaxy M31 
and the Milky Way Galaxy for wide regions from the cen¬ 
ters to the dark halo edges (figures 2 and 3). The GRC for 
the Milky Way was revised by adopting the most recent 
Solar rotation velocity and applying the same method as 
for M31. As stressed below, the dark halos of both galax¬ 
ies are well represented by the NFW density profiles. 
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Gomponent 

Parameter 

M31 

Milky Way 

Bulge 

Ob (kpc) 

1.35 ±0.02 

0.87 ±0.07 


Mb(10“MQ) 

0.35 ±0.004 

0.25 ±0.02 

Disk 

Od (kpc) 

5.28 ±0.25 

5.73±1.23 


Md(10“MQ) 

1.26 ±0.08 

1.12±0.40 

NFW Halo 

h (kpc) 

34.6±2.1 

10.7±2.9 


Po(10“^Mqpc“3) 

2.23 ±0.24 

18.2 ±7.4 


P8 kpc(10“^MQpC“3) 

6.36 ±0.70 

7.93 ±3.24 


= (in energy density: GeV cm“^) 

0.24 ±0.03 

0.30 ±0.12 


Mh:200(10“Mo) 

12.3 ±2.6 

5.7±5.1 


Mq) 

18.3 ±3.9 

7.3 ±6.7 

Total Mass 

fbftot : 200 ( 1 0 ^ ^ IHq ) 

13.9 ±2.6 

7.0 ±5.1 


fbftot : 385 ( 1 0 ^ ^ IHq ) 

19.9 ±3.9 

8.7±5.1 

Bulge 

Xt/N {Ri - i ?2 kpc) 

0.36 (0.0-20.0) 

3.5 (0.0-20.0) 

Disk 

xl/N {Ri - i ?2 kpc) 

0.33 (0.0-40.0) 

3.6 (0.0-40.0) 

Halo 

xl/N {Ri - i ?2 kpc) 

0.25 (0.0-385.0) 

3.0 (0.0-385.0) 


t .Mh: 200 : .^h: 385 j-^tot :200 ^nd Mtot :385 are dark halo and total masses within R = 200 and 385 kpc, respectively; pskpc 
is a local value at i? = 8 kpc both in mass and energy densities; Ri and i ?2 are start and end radii for fitting. 


By the least-y^ fitting to the obtained GRCs up to ra¬ 
dius of 385 kpc (figures 5 and 6), we determined the galac¬ 
tic parameters for the bulge, disk, and dark halo, as listed 
in table 2. Our result for M31 is consistent with those 
obtained by the other authors in the decade, as compared 
in table 3 and figure 9. 

4-2. The reality of NFW model 

The fitting result in figures 5 and 6 proves that the 
NFW profile (Navarro et al. 1995) can be a realistic ap¬ 
proximation to represent observed dark halos. We em¬ 
phasize that the analysis of GRGs covering regions as 
wide as several hundred kpc around the galaxies is es¬ 
sential to discriminate the right model among the three 
types of dark halo models: the isothermal model predict¬ 
ing flat rotation in the outermost region as V^ot oc const, 
the NFW model predicting slowly decreasing rotation as 
oc (In R/R)~^^^, and Plummer- or exponential-type mod¬ 
els predicting Keplerian decrease as oc at large radii. 

4-3. Similarity of the GRCs in MSI and the Milky Way 
and difference in the dark halo masses 

In figure 4 the GRC of M31 is compared with the re¬ 
vised GRG of the Milky Way, and the fitting results are 
compared in table 2. There is remarkable similarity of ro¬ 
tation curves inside the disk regions as well as in the dark 
halos in their shapes and amplitudes. However, the fitted 
dark halo mass of M31 is about twice that of the Galaxy 
and the parameters are also different, although the bulge 
and disk parameters are similar between the two galaxies. 

4.4- Revised GRC of the Milky Way 

The fitting accuracy in the Milky Way is not satis¬ 
factory, particularly for the bulge component (hgure 8), 
which is mainly due to the fact that the Galactic bulge 
is composed of multiple components, not well represented 
by the de Vaucouleurs law (Sofue 2013). However, we 


Table 3. Comparison of derived masses of dark halo and to¬ 
tal mass of M31 with other works in the decade. 


Authors 

Mass 

in IO^^Mq 

Evans et al. (2003) 

Mtot-.ioo 

7-10 

Ibata et al. (2004) 

Aftot:125 

7.5 ±1.3 

Geehan et al. (2006) 

A7h:20ot 

7.1 

Seigar et al. (2008) 

A7h:20ot 

7.3 

Lee et al. (2008) 

Mtot:100 

19.0±1.3 

Chemin et al. (2009) 

f^h:160 

10.0 

Corbelli et al. (2010) 

A7h:20ot 

13 ±3 

Watkins et al. (2010) 

Mtot:300 

14 ±4 

van der Marel et al. (’12) 

-^tot:Vir(385) 

17.2 ±2.5 

Tollerud et al. (2012) 

A7tot:139 

8±4 

Tamm et al. (2012) 

Arh:200t 

11.3-12.7 

Fardal et al. (2013) 

Mtot:200 

19 ±5 

Veljanovski et al. (2014) 

A7h:200t 

12-16 

This work (2015) 

A7h:200t 

12.3 ±2.6 

— 

A7h:385 

18.3 ±3.9 

— 

Aftot:200 

13.9 ±2.6 

— 

Aftot:385 

19.9 ±3.9 

Average < Mh: 2 oo 


10.8±3.0 

Average < Mtot:ioo -385 > 


13.7±5.2 


used the common model to M31 in order to compare the 
global structures of the two galaxies. Hence, the inner¬ 
most structure of the Milky Way may not be taken serious 
here, while it does not affect the result for the disk and 
dark halo. 

The here determined scale radii and masses for the 
Milky Way are systematically greater than those obtained 
in our previous work (Sofue 2013). This is because of the 
adopted larger rotation velocity of the Sun, which causes 
significant increase of the disk rotation velocities, as well 
as to the revised correction factor of radial velocities of 
non-coplanar objects. 

















6 


Yoshiaki Sofue 


[Vol. 


V 
(knVs) 


(a) 



(a) 



V 
(km/s) 


(b) 



(b) 



Fig. 5. Least-x^ fit of the GRC of M31 by the bulge, disk 
and dark halo components in (a) linear and (b) logarithmic 
scalings. 


Fig. 7. (a) Values of /N for M31 as functions of a^, <idi 
and h around the least values marked by circles. From left to 
right: bulge, disk, and dark halo, respectively, (b) Same, but 
for the masses Mb, M^ and Mh: 385 - 


V 
(km/s) 


(a) 



V 
(km/s) 


(b) 



Fig. 6. Same as figure 5, but for the revised new GRC of 
Milky Way using the same procedure as for M31. 


(a) 


(b) 



cn 

IE 

O 




10 ” 

Mass (Msun) 


10 '* 


Fig. 8. Same as figure 7, but for the Milky Way using the 
newly determined revised GRC. 
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R(kpc) 


Fig. 9. Enclosed dark (filled) and total masses (open cir¬ 
cles) within R of M31 obtained by the authors listed in table 
3. Upper (dashed) and lower (fill) lines show Mtot(R) and 
Mh(/?), respectively, calculated using the fitted parameters 
of M31. Big diamonds on the calculated lines indicate the 
present result. 

4-.5. Baryonic fraction 

The ratio of the total mass of the bulge and disk to the 
dark matter mass, T = (Mb + Md)/Mh: 385 , is a measure 
of the baryonic to dark mass ratio. The values in table 
2 yield the ratio of F ^ 0.087 for M31, and 0.25 for the 
Galaxy. These values may be compared with the cosmo¬ 
logical value of 0.19 (=4.5%/24%) from the WMAP ob¬ 
servations (Spergel et al. 2003). This implies that M31 is 
dark mater dominant compared to the cosmological value, 
whereas the Galaxy is baryon exceeded. 

4 . 6 . Relation to the Local Group 

The systemic velocity of M31 with respect to the Milky 
Way of ~ 100 km s“^(Courteau and van den Bergh 1999) 
can be translated to a statistically possible mutual ve¬ 
locity by multiplying to yield ^ 170 km s“^, if we 
assume that the two galaxies’ motions are random. In 
order for the two galaxies to be bound, a reduced total 
mass of the whole system is required to be greater than 
~ 5 X IO^^Mq. Since both M31 and the Galaxy are far less 
massive, ~ 3 x IO^^Mq in total, than this binding mass, 
there remain two possibilities either that the Local Group 
contains dark mass as massive as ^ 2 x 1 O^^M 0 , or that 
the two galaxy systems are not gravitaionally bound. 

In the former case, our view that the Local Group 
can be recognized as one system may not be changed. 
However, the dark mass must be concentrated around the 
center of mass of the Local Group, so that the M31 and 
the Galaxy groups, both associated with own dark halos 
and embedded satellite galaxies and globular clusters, are 
not tidally disrupted. 

In the latter idea, we need no further assumption of 
dark mass in the Local Group, but it is required that M31 
and the Galaxy and their satellite galaxies have been two 
individual bound systems since the formation. In fact, 
Sawa and Fujimoto (2005) proposed a model that dwarf 
galaxies of the Local Group are bound either to M31 or 
the Milky Way, composing two individual bound groups, 
and the two groups are orbiting around each other and 
tidally interacting. 
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Appendix 1. Pseudo rotation velocity in 
pressure-supported rotating system 

If there exists systematic rotation Vrot in the halo, the 
velocity H = |y | of a satellite galaxy and/or globular clus¬ 
ter (hereafter, particle) is expressed as 

y' = (yrot + yran)' = y,ot + 2yrot • yran + KL, (Al) 

where yran is the velocity corresponding to pressure- 
support term. When y^ is averaged around the galaxy, 
the crossing term disappears because of randomness of 
yran and axisymmetry of yrot, yielding 

< y2^ > ^ ^ (A2) 

The Hrst term is related to the apparent rotation velocity 
projected on the sky, Vz-.iot as 




Yoshiaki Sofue 


[Vol. , 

^j2:rot ~ ^ COS 0 . (A3) 

Here, i is the inclination angle of the rotation axis, and 9 
is the azimuthal angle of the particle from the major axis. 

Knowing that i and 9 are independent and < cos^9 >=1/3 
by averaging over 9 from 0 to 7r/2, we have 

< Kot >= < ^Lot > ■ (A4) 

Sin i 

Replacing < > by 3 < vl.^ >, where 'y^:ran is the 

z-directional component of Vran, we obtain 

. (A5) 

sm i 

Writing the observable z-directinal velocity as Vz = 

■Cziran ± t'zirot, and remembering that the crossing term 
< i2vz:TanVz:Tot > reduces to zero for randomness of sign 


of 'Czrran, we obtain 

< vl >=< > . (A6) 

Hence, we have 

< >= 3 <vl > +3cot^j < 'y^-rot > • (A7) 

For an edge-on case, i = 90°, we have 

<V^ >=3<vl>, (AS) 

the same as for random motion (Limber and Mathews 

1960). If we assume the same inclination as the main disk 
of M31, i = 77° to 78°,we have 

<V'^ >=3<vl> +a< -yf.rot > (A9) 

with a ~ 0.15. Since < > is smaller than < v'^ >, we 

may neglect the second term, obtaining 

<V^ >-3<vl> . (AlO) 


Finally, we define the pseudo rotation velocity as 


V = V<K2 > = y'S < >. 


(All) 
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